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Abstract 

In the minimal supersymmetric standard model (MSSM), a strongly first-order elec- 
troweak phase transition (EWPT) is only possible in a confined parameter region where one 
of the scalar top quarks is lighter than the top quark and the other one is as heavy as the 
SUSY breaking scale. If the MSSM is enlarged to accommodate vector-like quarks and their 
superpartners, we find that the strongly first-order EWPT is possible without requiring light 
scalar top quark at the one-loop level, in the limit where the lightest scalar Higgs boson of 
the MSSM behaves like the Higgs boson of the standard model and the other Higgs bosons 
are all as heavy as the SUSY breaking scale. The strength of the first-order EWPT is found 
to be dependent on the mass of the lightest neutral Higgs boson and the mixing effects of 
the vector-like scalar quarks. 
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I. Introduction 



Sakharov has suggested several decades ago that the baryon asymmetry of the universe may 
be dynamically generated if three conditions are met [1]. The three Sakharov conditions are 
the violation of the baryon number conservation, the violation of both C and CP, and the 
deviation from thermal equilibrium. Among them, the deviation from thermal equilibrium is 
the most intricate and challenging condition for theoretical models to satisfy. In order to ensure 
sufficient deviation from thermal equilibrium, the electroweak phase transition (EWPT) should 
be strongly first order, since otherwise the baryon asymmetry generated during the EWPT 
would subsequently disappear. 

In general, the electroweak symmetry of the universe is spontaneously broken as the universe 
is cooled down from very high temperature to the present temperature. In between the cooling 
process, the shape of the potential of the universe at a critical temperature would have two local 
minima for the first-order EWPT. One minimum corresponds to the original symmetric-phase 
state and the other one the broken-phase state. The strength of the EWPT is measured by 
comparing the vacuum expectation value (VEV) of the Higgs field at the broken-phase state 
with the critical temperature. If the critical VEV is larger than the critical temperature, the 
first-order EWPT is said to be strong. 

A number of supersymmetric models have been investigated within the context of the EWPT 
for the electroweak baryogenesis. Several authors have investigated the MSSM to show that a 
strongly first-order EWPT is only possible in a confined parameter region where a scalar top 
quark is lighter than the top quark and the other scalar top quark is as heavy as the SUSY 
breaking scale of 2 TeV [2-5]. Then, in order to avoid the light scalar top quark, other models 
extended from the MSSM have also been investigated for the possibility of a strongly first-order 
EWPT, such as the non-MSSMs [6-19] and the MSSM with four generations of fermions [20,21]. 

The non-MSSMs are different from the MSSM in the sense that there is a trilinear term of 
the singlet Higgs field accounting for the effective /i term. The trilinear term between Higgs 
singlet and doublets of the non-MSSMs allows the mass of the scalar top quark to be sufficiently 
large. Thus, the non-MSSMs have some parameter regions where a sufficiently strong first-order 
EWPT may occur for electroweak baryogenesis, without requiring a light scalar top quark. 
Recently, it has been found that if the fourth generation of fermions is introduced to the MSSM, 
it is possible to have a strongly first-order EWPT without a light scalar top quark in terms of 
the thermal contribution of the scalar quarks of the fourth generation [20,21]. 

In this article, we study another model, the MSSM with vector-like quarks within the context 
of the electroweak baryogenesis. It has been proposed some years ago [22,23], and its Higgs sector 
has been extensively investigated recently [24-28]. There are three ways to introduce the chiral 
superfields to the MSSM, maintaining the perturbative gauge unification with the masses of new 
extra chiral supermultiplets at the TeV scale [28]. Among them, we pay attention to the specific 
model, where two vector-like quarks are present. This model is similar to the QUE model in 
Ref. [28], except that we do not consider an additional chiral lepton superfield. 

Our model is therefore actually the MSSM with two additional chiral superfields, Q and 
which are defined to transform under S'C/(3)c x SU(2)l x U{1)y as 

Q = (3,2,1/6), Q= (3,2,-1/6) , 

= (3,1,2/3), = (3,1,-2/3) . (1) 

These additional chiral superfields introduce new physical vector-like quarks, namely, two charge 
-|-2/3 quarks Ti, T2, and one charge —1/3 quark B, together with their scalar superpartners, Tj 
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{i = 1,2,3,4) and B {i = 1,2). 

The presence of the vector-hke sector in the MSSM alters the phenomenology of the MSSM 
with respect to the EWPT. Our study shows that in our model a strongly first-order EWPT is 
possible at the one-loop level, without requiring a light scalar top quark. 

II. Higgs potential 

The superpotential of our model may be written as 

W = Wmssm + MqQQ + MuUU + YunuQU -YaUdQU , (2) 

where Wmssm is the superpotential of the MSSM, Mq and Mfj are the masses of the vector-like 
quarks, and and Y^ are the Yukawa coupling coefficients to the weak hypercharge —1/2 and 
1/2 Higgs doublets HJ = (i?)^,F^) and = {H+,H^) in the Higgs superfields Tiu and Hd, 
respectively. 

The tree-level Higgs potential of the MSSM is given as 

= ml\H^\^ + mj\Hd\^ - [ml^H^Hd + H.c.) + ^(^i' + ^i) {Wu? - l^^dP)' , (3) 

where rriu, m^, and mud are the mass parameters, gi and §2 are the gauge coupling coefficients 
for U{1) and SU{2), respectively. Two of the mass parameters, nid and m^, may be eliminated 
using the minimum equations with respect to the neutral Higgs fields. 

After electroweak symmetry breaking, as the neutral Higgs fields develop VEVs, Vd = {Hd) 
and Vu = {Hu), five physical Higgs bosons emerge: two neutral scalar Higgs bosons, h and H 
{rrih < mn), one neutral pseudoscalar Higgs boson, A, and a pair of charged Higgs bosons, H^. 
We assume that the VEVs as well as the parameters of the Higgs potential are all real. Note 
that the tree-level mass of h can be expressed in terms of just two parameters, rriud and tan /?, 
where tan /3 = Vu/vd- The tree-level mass of h is smaller than the Z boson mass. 

Let us consider the decoupling limit where h remains light at the electroweak scale and all 
the other Higgs bosons become very heavy with masses of the order of the SUSY breaking scale 
or a few TeV. In this limit, the Higgs sector reduces to the Higgs sector of the Standard Model 
(SM) at the electroweak scale and the coupling coefficients of h to gauge bosons and fermions 
are identical to the couplings of the SM Higgs boson, (p. Thus, in the decoupling limit, one 
cannot distinguish the phenomenology of h from that of 0, hence the same experimental lower 
bound on the mass of h as the SM Higgs boson which is about 114.4 GeV as set by the LEP2 
experiments [29]. 

Hereafter, we confine ourselves in the decoupling limit and denote the lighter neutral scalar 
Higgs boson as </>: 

</> = cos /SReiH^) + sin PRe{H^) . (4) 
In terms of (j), the tree-level Higgs potential at zero temperature may be written as 

W,O) = -^0^ + ^<AS (5) 

where tuq is related to the mass of the lighter neutral scalar Higgs boson at the tree level and 
= idi +52) cos^ 2/3/4 represents the quartic Higgs self-coupling coefficient. 
Let us first study the Higgs potential of our model at zero temperature. The tree-level 
Higgs potential of our model at zero temperature is essentially identical to the tree-level Higgs 
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potential of the MSSM. However, the one-loop level Higgs potential of our model differs from 
the corresponding MSSM Higgs potential, since additional contributions from the loops of the 
vector-like quarks and the vector-like scalar quarks enter into the radiative corrections to the 
Higgs potential of our model. At zero temperature, the one-loop effective potential of our model 
may be written as 



y(0,o) = 1/0(0,0) + yi(<^,o) , 



(6) 



where Vi(i;^>, 0) is the radiative corrections at zero temperature. It is given by the effective 
potential method as [30] 



nimf{(l)) 
647r2 



log 



A2 



(7) 



where I stands for the participating particles, namely, the electroweak gauge bosons W and Z, 
top quark t and scalar top quarks ti and t2, as well as the vector- like quarks Ti and T2 and the 
vector-like scalar quarks Tj = 1, 2, 3, 4), A is the renormalization scale of the modified minimal 
subtraction scheme, set as mz, and n/ is the degree of freedom for the participating particles: 
nw = nz = 3, nt = —12, n^^ = 6 (z = 1, 2), n^i = n-Ta = —12, rif = 6 (i = 1, 2, 3, 4). 

Now, we calculate the tree-level masses of the relevant particles which will be used to obtain 
the one-loop effective potential. The masses of W boson, Z boson, and top quark, respectively, 
at the tree-level are given as 



2 / 2 2\ 

^ww) = ' ^zw) = — 2 — ^ 



hf sin ^2^2 



(8) 



We assume that the lighter scalar top quark is predominantly right-handed. As is well known, 
the scalar top quarks contribute to the tree-level Higgs sector considerably. The masses of the 
scalar top quarks are given as [2,4], 



m?^ (0) = + 



D 



1 



1 + 



m 
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where mg^ and mj/g are the SUSY breaking mass parameters, -D^ ((/>) and Dj^{(l)) are the D-term 
contributions arising from the SM gauge bosons, defined as 



DlicP) = ^^:^cos2P<l>\ 2^1(0) = I cos --^2 



(10) 



and At = At — ^cot/3, where At is the trilinear SUSY breaking mass parameter. Note that the 
difference in size between D'j^{(p) and D'^^cf)) causes a relatively small mass splitting between 
ti and ti, whereas the opposite sign of At in the above expressions indicates that the mass 
eigenstates of scalar top quarks are mixtures of the weak eigenstates. 

We also calculate the masses of the vector-like quarks and the vector-like scalar quarks at 
the tree level. The mass matrix for the vector-like quarks at the tree level is given as M^Mj, 
[28], with 



Mt((/)) 



Mq 
Yd cos / 



Mu 



:ii) 
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The eigenvalues of the mass matrix are the squared masses of the vector-hke quarks. If we take 
Mq = Mjj for simphcity, they are given to a good approximation as: 



^Q-^ ^us fJfp ■ (12) 
The mass matrix for the vector-hke scalar quarks are expressed as a 4 x 4 matrix, 



/ ml + Dl{<t>) 


V Yu sin f3(j)Au 







Yd cos P(j)Aa 



m 



Q 



Yd cos (3 <f> Ad 



Yu sin 13 (f>Au \ 



Di 







,(13) 



where bq and bu are the soft SUSY breaking mass, iuq, itiq, mu and m^; are the soft SUSY 



breaking masses, A^ = A^/Y^ — ficot (3 and Ad = Ad/Yd 
trilinear soft SUSY breaking mass parameters and 



Mi 



MtM. 







M\,Mt 



fj. tan /3, with A^ and Ad being the 
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In subsequent calculations, we set bq = bjj = 0, Mq = Mu, and adapt some approximations for 
the mass matrices, following Ref. [28]. 

The squared masses of the four vector-like scalar quarks are given as the eigenvalues of their 
mass matrix. We would like to reduce the number of free parameters, by absorbing Mq (= Mjj) 



mf 



into mq, uiq, mu and rriu- Thus, we may redefine MQ + mjy as m^, Mq + uiq as m^, Mq + 
as nifj, and Mq + tuq as m^. Then, the squared masses of the vector-like scalar quarks are 
given in simplified expressions as 



(0) 



m|^((/)) 



ml + 



R\ 



+ sin^ , 



2 
Q 



m'Q + Dii<P)+Y^'sinU 
+ Y? cos^ , 



1 + 



m 



mfj 
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m% 



m 



Q 



+ Yi cos2 , 



1 + 



i3 



my 



(15) 



We note that we keep the D-term contributions, both in the scalar top quark masses and in 
the vector-like scalar quark masses. Due to the presence of Di and Dr, these scalar quarks are 
not degenerate in mass in general. If both Du = Dji = and Ag = (q = t, u, d), the scalar top 
quarks as well as the vector-like scalar quarks would be degenerate in mass. 

We would like to note that the strength of the electroweak phase transition tends to decrease 
when the mixings between scalar quarks are taken into account. Thus, for the first order EWPT 
to be strong, the mixings among the scalar quarks are not favored. Nevertheless, we assume that 
Aq {q = t, u, d) are not zero, implying that there are mixings among scalar quarks. Later, we will 
discuss on the effect of the mixings among the scalar quarks on the strength of the first-order 
EWPT. 
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The renormalized parameter ttiq in the Higgs potential can be expressed in terms of other 
parameters. Prom the minimum condition of the first derivative of V{4>, 0) with respect to (j), 
niQ at the one-loop level is expressed as 
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1 



cos^ 2/3 
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(16) 



where the first term on the right-hand side is obtained from Vo(</), 0) and the remaining other 

terms are obtained from Vi{(l),0). Also, v = \Jv'^ + = 175 GeV and the field-dependant 
tree-level masses for the relevant particles are evaluated at the electroweak symmetry breaking 
scale. The mass of (p the one-loop level is obtained by calculating the second derivative of 
y(</!), 0) with respect to (j), evaluated at (p = v, as 



rrij 



m|cos2 2/3 + |5lllog 



3m^ 




3m^ 



1 I !I!± 
47r2?;2 I A2 , 



(17) 
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where the first term on the right-hand side is the tree-level value and the remaining terms are 
the radiative corrections. The field-dependant tree-level masses for the relevant particles are 
evaluated at the electroweak symmetry breaking scale, v. 

Let us now study the Higgs potential of our model at finite temperature. At temperature T, 
the full effective Higgs potential at the one-loop level may be written as 

y(0,r) = yo(0,o) + yi(</.,o) + yi((A,T) , (is) 

where Vi{(f),T) is the the radiative correction at finite temperature. Employing the effective 
potential method [31], we have 

(19) 



V^{(t>. T) = Y^'^ dx x^ log 1 ± exp (^-\/x2 + m2(0)/r2 



where I = W, Z, ii, {2, Ti, T2, % {i = 1,2,3,4), and the negative sign is for bosons and the 
positive sign for fermions. Note that the thermal effects due to W boson, Z boson, top quark, 
vector- like quarks and their scalar partners are taken into account. The possibly degenerate 
global minimum of V((/), T) is defined as the vacuum at T. The explicit expression for Vi{(f),T) 
is obtained by the exact numerical integration, instead of the high-temperature approximation. 



III. Numerical Analysis 

It is easy to confirm that the above Higgs potential at the one-loop level becomes symmetric 
at very high temperature. Also, one may check that the electroweak symmetry of the Higgs 
potential is broken at zero temperature. In between the two extreme temperatures, there is an 
intermediate temperature, called the critical temperature Tc, at which the above Higgs potential 
has a doubly degenerate minimum. Thus, at Tc, one minimum occurs at = 0, which is the 
trivial minimum, and the other minimum takes place at (j) = Vc, the critical VEV. The trivial 
minimum represents the vacuum of the symmetric phase of the universe and the minimum of 
<j) = Vc represents the vacuum of the broken phase. The EWPT may take place from the 
symmetric phase vacuum to the broken phase vacuum. 

The critical temperature as well as the critical VEV is important for the successful EWPT, 
since they determine the strength of the first-order EWPT. Indeed, Vc defined in this way might 
be a reasonable characterization of the strength of the phase transition. The strength of the 
first-order EWPT, which is driven by a non-equilibrium (out-of-equilibrium) condition, depends 
on the relative size between the critical temperature Tc and the critical VEV Vc- In general, the 
first-order EWPT is said to be strong if Vc> Tc- 

For our numerical analysis, we fix the values of some parameters. We set A = mz for the 
renormalization scale and mj = 175 GeV for top quark mass. Also, we take fng^ = 6 x 10^ GeV^, 
= 10^ GeV^, Ji^ = 6 x 10^ GeV^, which appear in the scalar top quark masses, Mq = 700^ 
GeV2, m}j = ml = 200^ GeV^, m| = m| = 2.4 x 10^ GeV^, and Al = Al = 2.16 x 10^ GeV^, 
which appear in the masses of the vector-like scalar quarks. 

The Yukawa coupling coefficients of the vector-like quarks in our model should satisfy the 
renormalization group (RG) equations in order to avoid possible Landau poles at high energy 
scale, as discussed by Moroi and Okada [22,23]. The RG equations for and Yd are solved 
numerically by using the Runge-Kutta method, where no Landau pole is present. The results 
for their evolutions are shown in Figs, la and lb, where, for simplicity, we take the same initial 
values for Y^ and Y^ at the electroweak scale: Yu = Y(i = 0.87, 0.9, 1.0, 1.1, and 1.13. 
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Notice that if the initial values for and Yd are not larger than 1.0, both Yu and Yd become 
smaller for higher energy scale. However, if they start from initial values larger than 1.0, as 
Figs, la and lb show, they tend to increase as the energy scale increases. This is because their 
renormalization beta functions are dominated by the positive contributions from themselves 
which are larger than the negative contributions from the gauge couplings. 

With these values, we calculate the value of V{(j), T) as a function of (p at given temperature 
T, such that V{0,Tc) = V{vc,Tc) = and 

dV _9y 

These two conditions establish the doubly degenerate vacua of the symmetric phase and the 
broken phase. 

In this way, we obtain Fig. 2 for tan/3 = 5 and 1^ = 1^ = 1, at r = 157.58 GeV, where 
V{(p, T) shows the typical shape for the first-order EWPT. One can see that there are two 
minima for V{4>, T), at around (j) = and (j) = 258 GeV. Since the values of V{(p, T) at these two 
minima are equal, they are doubly degenerate vacua with a potential barrier between them. The 
first-order EWPT may take place from one vacuum to the other one through thermal tunneling. 
Therefore, (f) = 258 GeV is the critical VEV for the vacuum of the broken phase and T = 157.58 
GeV is the critical temperature Tc that allows the phase transition from the symmetric-phase 
state to the broken-phase state. From these values, we obtain Vc/T^ ~ 1.637, which certainly 
implies that the first-order EWPT shown in Fig. 2 is strong. Consequently, Fig. 2 suggests that 
there is at least one set of parameter values of our model at the one-loop level which allows the 
possibility of a strongly first-order EWPT. 

For the set of parameter values of Fig. 2, the masses of participating particles are evaluated. 
The mass of the lightest scalar Higgs boson is calculated to be = 151.41 GeV, whereas the 
masses of the other Higgs bosons are all very heavy in the decoupling limit. The masses of the 
scalar top quarks are obtained as m^-^ = 190.67 GeV and m^^ = 794.35 GeV. Note that both of 
two scalar top quarks are heavier than top quark. In other words, our model does not require a 
light scalar top quark in order to realize a strongly first-order EWPT. 

Vector-like quarks and scalar quarks are characteristic members of our model. Their masses 
are calculated to be ttt-Ti = 720.72 GeV and = 700.84 GeV. The masses of the vector- 
like scalar quarks are mrif^ = 200.62 GeV, m^,^ = 1566.76, m^,^ = 203.28, and ruf^ = 1550.97. 
These masses might affect the precision electroweak observables, due to the virtual corrections 
at the one-loop level. As noticed by Ref. [28], the most important contributions from the 
vector-like quark sector to the precision electroweak observables may be described in terms of 
two parameters, S and T, which have been proposed by Peskin and Takeuchi [32]. We follow 
the method of Ref. [28] to calculate the contributions from the vector-like quark sector to the 
electroweak gauge boson self energies. We find that the corrections are practically negligible: 
A5 = 0.00187 and AT = 0.00686 for the parameter values in Fig. 2. Thus, the masses of the 
vector-like quark sector obtained from the parameter values in Fig. 2 are consistent with the 
electroweak precision test. 

Now, we would like to show that the above set of parameter values is not singular. We 
examine other regions of the parameter space to find that a large number of sets of parameter 
values within the parameter space allow the strongly first-order EWPT. As an illustration, 
we vary the value of tan/3, while fixing the values of other parameters, and repeat the above 
procedure. Our result is listed in Table 1, where one may note that a wide range of tan/3 indeed 
allow the desired phase transition. 
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Table 1: The critical temperature and the critical VEV for the strongly first-order EWPT for 
some values of 1 < tan/3 < 40. The values of the other parameters are the same as Fig. 2, in 
particular, Al = = 2.16 x 10^ GeV^ and Y^ = Yd = 1. 



tan /3 


1 


1.5 


2 


5 


10 


30 


40 


Tr (GeV) 


146.19 


150.83 


153.34 


157.58 


158.43 


158.69 


158.71 


Vc (GeV) 


289 


285 


280 


258 


255 


254 


253 


(GeV) 


107.84 


116.27 


127.81 


151.41 


156.15 


157.64 


157.73 


Vc/Tc 


1.9768 


1.8895 


1.8260 


1.6372 


1.6095 


1.6006 


1.5941 


(GeV) 


193.81 


192.51 


191.77 


190.67 


190.48 


190.41 


190.41 


(GeV) 


796.04 


795.33 


794.94 


794.35 


794.24 


794.21 


794.21 


(GeV) 


710.85 


714.98 


717.28 


720.72 


721.33 


721.51 


721.53 


mra (GeV) 


700.85 


706.69 


704.36 


700.84 


700.21 


700.02 


700.01 


ruf^ (GeV) 


203.79 


202.47 


201.73 


200.62 


200.42 


200.36 


200.36 


rrif^ (GeV) 


1558.55 


1561.97 


1563.89 


1566.76 


1567.26 


1567.42 


1567.43 


rrif,^ (GeV) 


203.79 


203.57 


203.45 


203.28 


203.24 


203.23 


203.23 


mf, (GeV) 


1558.55 


1555.40 


1553.63 


1550.97 


1550.50 


1550.36 


1550.35 



After establishing the possibility of the strongly first-order EWPT, we are now interested in 
the behavior of the strength of the first-order EWPT with respect to the mass of the lightest 
scalar Higgs boson. It has been observed in both the SM and the MSSM that the strength of 
the first-order EWPT always increases as the lightest scalar Higgs boson mass decreases. Also, 
in the non-MSSMs with a Higgs singlet, a stronger first-order EWPT is allowed for a lighter ^i, 
where 5i is the lightest scalar Higgs boson [6-19]. Thus, one can expect that in our model too 
the strength of the EWPT is enhanced as the lightest Higgs boson becomes light. 

In Table 1, one can note that increases as tan (3 increases from 1 to 40. At the same time, 
Vc/Tc decreases as tan/3 increases from 1 to 40. Thus, for the set of parameter values in Table 
1, one may notice the tendency of stronger phase transition for smaller m^. In other words, 
the behavior of our model is consistent those of the SM, the MSSM, or the non-MSSMs in the 
sense that the smaller mass of the lightest Higgs boson may indeed enhance the strength of the 
first-order EWPT. 

Another point we would like to discuss is the role of the Higgs cubic term that appears in the 
power expansion of the one-loop Higgs potential at finite temperature, in the high-temperature 
approximation. The size of this Higgs cubic term is essential to the strength of the first-order 
EWPT. The coefficient of this term depends on a number of bosons such as the SM gauge 
bosons, the scalar top quarks, and the vector-like scalar quarks. The SM gauge bosons may 
take part in the coefficient of the (jy^ term in the SM. However, in the MSSM, its contribution 
to the coefficient of the (ji^ term is not so much important as the contributions from the scalar 
top quarks, as long as the scalar top quarks are sufficiently light. 

If the masses of the scalar top quarks are very large (1 or 2 TeV) in the MSSM, its contribution 
to the strength of the first-order EWPT through the Higgs cubic term would be negligibly small. 
On the other hand, if they are sufficiently light at the electroweak scale (for example, about 165 
GeV between the top quark mass and the experimental lower bound of about 155 GeV), it might 
significantly contribute to the strength of the first-order EWPT. Thus, in the MSSM, at least 
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one of the scalar top quarks need to be lighter than top quark in order to achieve a strongly 
first-order EWPT. 

The strength of the first-order EWPT in our model depends also on the vector-like scalar 
quarks. As shown in Eq. 15, the pattern of the masses of the vector- like scalar quarks may be 
represented by a sum of the soft SUSY mass and the Yukawa term, up to a constant. The 
terms are constant if the gauge coupling coefficients remain as constant. 

If the Yukawa coupling coefficients are small (< 0.1) and the soft SUSY masses are large 
(> 1 TeV), the masses of the vector-like scalar quarks would be large (> 1 TeV), and the 
loop contributions of the vector-like scalar quarks to the Higgs cubic term would be very small. 
Thus, in this case, the vector-like scalar quarks suppress the strength of the first-order EWPT, 
by suppressing the coefficient of the cj)^ term. In fact, in the symmetric phase state where (p 
becomes zero, the vector-like scalar quarks do not become massless but remain heavy at TeV 
scale. In other words, large soft SUSY masses prevent the appearance of the Higgs cubic term 
in the expansion of the effective thermal potential. 

On the other hand, if the Yukawa coupling coefficients are large (~ 1) and the soft SUSY 
masses are small 200 GeV), some of the vector- like scalar quarks would be light (~ 200 
GeV), In this case, the loop contributions of the vector-like scalar quarks to the Higgs cubic 
term would be large, and thus enhance the strength of the first-order EWPT. Actually, our 
choice of parameter values are large Yukawa coupling coefficients and small soft SUSY masses 
such that a strongly first-order EWPT is possible. 

The situation is quite similar in our model, where there are additional scalar fields, namely, 
vector-like scalar quarks. When Vi{<j),T) in our model is expanded in powers of (j) in the high- 
temperature approximation, there is a term proportional to c/)^, arising from the vector- like scalar 
quarks as well as from the SM gauge bosons and scalar top quarks. The additional contribution 
due to loops of vector-like scalar quarks certainly enhance the strength of the first-order EWPT, 
thus relieving the necessity of a light scalar top quark in our model. 

Next, let us study the effect of the mixing in the scalar quarks on the strength of the first-order 
EWPT. The mass eigenstates of scalar quarks are given as mixtures of the weak eigenstates, 
since the mass matrix for them is not diagonal in general. The size of the mixing between them 
is determined by the size of the off-diagonal matrix elements. In our model, the mixing between 
the scalar top quarks is determined by Af and the mixing among the vector-like scalar quarks 
is determined by = A^. We show in Table 2 the result of varying these parameters, while 
keeping other parameters fixed. 

The result in Table 2 suggests that the first-order EWPT becomes much stronger if the value 
of the parameter for the mixing among vector-like scalar quarks is small. 

IV. Conclusions 

By introducing the vector-like quarks and scalar quarks to the MSSM, we study the possibility of 
a strongly first-order EWPT, at the one-loop level. The Higgs sector of our model is essentially 
identical to the MSSM Higgs sector, while there are additionally two vector-like quarks and four 
vector-like scalar quarks in our model. At the one-loop level, we take into account the thermal 
effects due to loops of W boson, Z boson, top quark, scalar top quarks, vector-like quarks, and 
vector-like scalar quarks. 

We find that a number of parameter sets in the reasonably wide parameter space of our 
model allow the strongly first-order EWPT, thus satisfying one of the Sakharov conditions 
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Table 2: The critical temperature and the critical VEV for the strongly first-order EWPT for 
two different values of Af^ = A\. The values of the other parameters are the same as Table 1, in 
particular, tan/3 = 5. Thus, the first column is identical to the fourth column of Table 1. 



Ai = Ai (GeV^) 


2.16 X 10^ 


1.2 X 10^ 


Tc (GeV) 


157.58 


177.65 


Vc (GeV) 


258 


522 


(GeV) 


151.41 


122.84 




1.6372 


2.9383 


m-u (GeV) 


190.67 


190.67 


m-t, (GeV) 


794.35 


794.35 


rriTi (GeV) 


720.72 


720.72 


(GeV) 


700.84 


700.84 


nif^ (GeV) 


200.62 


228.09 


mf^ (GeV) 


1566.76 


1562.99 


nif^ (GeV) 


203.28 


203.28 


(GeV) 


1550.97 


1550.82 



for electroweak baryogenesis. The strength of the first-order EWPT is very important for the 
successful electroweak baryogenesis. In the MSSM, a strongly first-order EWPT is guaranteed 
by a light scalar top quark; whose mass should be smaller than the top quark mass. In other 
words, it is experimentally challenging for the MSSM to discover a scalar top quark whose mass 
is smaller than the top quark mass, if the first-order EWPT should be strong enough. On the 
other hand, in our model, the masses of the scalar top quarks need not be smaller than top 
quark. 

We work in the limit where all the Higgs bosons except for the lightest scalar Higgs boson are 
as heavy as the SUSY scale of a few TeV. In this limit, the lightest scalar Higgs boson behaves just 
like the SM Higgs boson whereas the other Higgs bosons are decoupled from the phenomenology 
at the electroweak scale. We note that, in the MSSM, the strength of the first-order EWPT 
increases as the mass of the pseudoscalar Higgs boson increases. Thus, the decoupling limit may 
be considered as an optimal situation for the first-order EWPT to be sufficiently strong, since 
the mass of the pseudoscalar Higgs boson is already very large. However, the decoupling limit 
alone may not be enough for the MSSM to achieve the desired phase transition, since in this 
limit the MSSM Higgs sector resembles the SM Higgs sector and it has already been observed 
that the SM Higgs sector cannot successfully generate the baryon asymmetry of the universe. 
Hence, the MSSM requires a different mechanism in addition to the decoupling limit for the 
strongly first-order EWPT. Our study suggests that the introduction of vector-like quark sector 
may be a possible mechanism. The existence of the vector-like quarks and scalar quarks not only 
enhance the strength of the first-order EWPT in our model, but also alleviate the requirement 
that the mass of a scalar top quark should be lighter than the top quark mass. 

The strength of the first-order EWPT in our model is affected by the values of the relevant 
parameters. Our numerical analysis shows that it is enhanced as the mass of the lightest scalar 
Higgs boson decreases. Also, if we reduce the value of = A^, the parameters that account 
for the mixing among the vector-like scalar quarks, the strength of the EWPT increases. 
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Figure la: The RG evolutions for by using Runge-Kutta method for different initial values 
of Yu and at the electroweak scale. 
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Figure lb: The RG evolutions for by using Runge-Kutta method for different initial values 
of Yu and Y^ at the electroweak scale. 
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Figure 2: The values of the relevant parameters are mq,^ = 6 x 10^ GeV^, ?n^^ = 10^ GeV^, 
= 6x 10^ GeV2, = 700^ GeV^, mfj = mfj = 200^ GeV^, = m| = 2.4 x 10^ GeV^, 
and Al = Al = 2.16 x 10*^ GeV^, and Yu = Yd = 1.0 at the electrowak scale. 
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